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Abstract 

Topology is an important degree of freedom in characterizing electronic systems. Recently, 
it also brings new theoretical frontiers and many potential applications in photonics. However, 
the verification of the topological nature is highly nontrivial in photonic systems as there is no 
direct analog of quantized Hall conductance for bosonic photons. Here we propose a scheme of 
investigating topological photonics in superconducting quantum circuits by a simple parametric 
coupling method, the flexibility of which can lead to the effective in situ tunable artificial gauge 
field for photons on a square lattice. We further study the detection of the topological phases 
of the photons. Our idea employs the exotic properties of the edge state modes which result in 
novel steady states of the lattice under the driving-dissipation competition. Through the pumping 
and the photon-number measurements of merely few sites, not only the spatial and the spectral 
characters, but also the momentums and even the integer topological quantum numbers with 
arbitrary values of the edge state modes can be directly probed, which reveal unambiguously the 
topological nature of photons on the lattice. 
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INTRODUCTION 


Charged particles in two dimension exhibit integer quantum Hall effect (IQHE) when 
exposed to a perpendicular magnetic field [T], characterized by the quantized transverse 
conductances in transport experiments. This novel effect can be explained by the integer 
topological Chern numbers describing the global behavior of the energy bands Such 

topological insulating IQHE phase is robust against disorder and defects because the band 
topology remains invariant as long as the band gaps are preserved. Therefore, in the con¬ 
nection between the topologically nontrivial material and the trivial vacuum, there exist 
unavoidably the edge state modes (ESMs) spatially confining at the boundary and spec¬ 
trally traversing the band gaps [3]. The presence of these gapless ESMs thus serves as an 
unambiguous signature of the topological non-triviality of the bulk band structure. 

Recently, the concept of topology has been extended to circuit quantum electrodynamics 
(QED) lattice [5J 6J, where the electrons are replaced by microwave photons hopping between 
superconducting transmissionline resonators (TLRs) [7, |8j. While the idea of topological 
photonics was firstly developed in photonic crystals [9HI3], circuit QED enjoys the time- 
resolved engineering of a large-scale lattice at the single-site level mm- The demonstrated 
strong coupling between superconducting qubits and TLRs 0IH! further allows the effective 
photon-photon interaction P7)HT§| which can hardly be achieved in other physical systems, 
indicating prospective future of investigating strongly correlated photonic liquids [20, ;21j. 
As photons are charge neutral, there have been several proposals of synthesizing artificial 
magnetic fields on a TLR lattice, with predicted strengths much stronger than those in 
conventional electronic materials [221 - 125] . Nevertheless, the synthetic Abelian gauge field 
has not been implemented so far despite the extensive theoretical studies, partially due to 
the complicated circuit elements required in these schemes. In addition, the detection of 
the integer topological invariants has also been addressed in recent research [2BH2B], which 
is nontrivial in the sense that the Hall conductance measurement cannot be transferred to 
circuit QED due to the absence of fermionic statistics. 

Here, we propose a theoretical scheme of implementing topological photonics in a two- 
dimensional circuit QED lattice. The distinct merit of our proposal is that we couple the 
TLRs by parametric frequency conversion (PFC) method which is simple in experimental 
setup and feasible with state-of-the-art technology [2HHII2] • The lattice in our scheme is 
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formed by TLRs connected to the ground through the superconducting interference devices 
(SQUIDs) [29, ME El EH, where the tunable photon hopping with nontrivial phases between 
TLRs can be induced through the dynamic modulation of the SQUIDs, allowing the arbitrary 
synthesization of time- and site-resolved gauge fields on the square lattice (1X1183] . Moreover, 
with the driving-dissipation mechanism being employed, various quantities of the ESMs can 
be measured by the pumping and the steady-state photon number (SSPN) detection of only 
few sites on the lattice [3J2J. In particular, the integer topological winding numbers of the 
ESMs with arbitrary values can be directly probed through the realization of the adiabatic 
pumping process HEX]. Such measurement is equivalent to the measurement of the Chern 
numbers of the bulk bands and thus clearly examine the topological nontriviality of the 
photons. Furthermore, our detailed discussions show unambiguously that our proposal is 
very robust against various potential imperfection sources in experimental realizations due 
to the topological nature of the ESMs, pinpointing the feasibility with current level of 
technology. Being flexible for the extension to more complicated lattice configurations and 
the incorporation of effective photon correlation, our scheme serves therefore as a promising 
and versatile platform for the future investigation of various photonic quantum Hall effects. 


RESULTS 

The lattice. —We start with a square lattice consisting of TLRs with four different 
lengths placed in an interlaced form, as shown in Fig. [IJa). At their ends, the TLRs are 
commonly grounded by SQUIDs with effective inductances much smaller than those of the 
TLRs PHI] . Due to their very small inductances, the grounding SQUIDs impose low- 
voltage shortcuts at the ends of the TLRs. Therefore the lowest eigenmodes of the lattice 
can be approximated by the A/2 modes of the TLRs with their ends being the nodes, and 
the whole lattice can be described by 

T~Lg ^ ] uj r ala r . ( 1 ) 

r 

with a\/a r being the creation/annilhilation operators of the rth photonic mode and uj v being 
the eigenfrequency. We further specify the eigenfrequencies of the four kinds of TLRs as 
yellow-wo, blue-t^o + A, green-cuo + 3A, and red-wo + 4A, respectively, with(Uo/27r e [10,20] 
GHz and A/2-7T G [1,3/2] GHz. Such configuration is for the following application of the 
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dynamic modulation method and can be achieved through the length selection of the TLRs 
in the millimeter range [3DH321 EE] • 

We then consider how to implement on the TLR lattice the effective tight-binding Hamil¬ 
tonian 

Ht = T ^ al,a r e~ l9r ' r + h.c., (2) 

<r,r'> 

in the rotating frame of Hs- Here T is the uniform hopping amplitude, and 

#r'r = J A(x) • dx, (3) 

is the r —» r' hopping phase manifesting the presence of a vector potential A(x) through 
Peierls substitution [3j. For each plaquette of the lattice, the summation of the hopping 
phases around its loop has the physical meaning 


^plaquette y A(x) • dx 

i. e. the synthetic local magnetic field for the microwave photons. 

However, it is nontrivial to have complex hopping constants between TLRs because the 
physical coupling between two TLRs takes real coupling constants, no matter capacitive 
dug or inductive [3HI32]. We then consider the dynamic modulation method studied 
in recent experiments [3TH - ET2] . The grounding SQUIDs can be modeled as flux-tunable 
inductances and it is now experimentally possible to modulate the SQUIDs by a. c. magnetic 
flux oscillating at very high frequencies (the experiment-achieved range is typically 8 ~ 10 
GHz [291, 3BJ which is much higher than the following-proposed 1 ~ 6 GHz). Such a. c. 
modulation introduces a small a. c. coupling 


// B(x) 


dS, 


(4) 


%ac = ^ 7rf{t)(a r + <4)(a r / + al), (5) 

<r,r') 

in addition to the d. c. contribution of the SQUIDs which is irrelevant because the TLRs are 
largely detuned [35j. We then assume that the a. c. modulation of the grounding SQUIDs 
contain three tones with frequencies being A, 2A, and 4A. By bridging the frequency 
differences between the TLRs, the 2A/4A tones induce the vertical blue <=t green/red <=t 
yellow parametric hoppings, and the A tone establishes the horizontal yellow blue and 
green red parametric hoppings [3J5J. When experiencing the PFC process between the 
TLR sites, the microwave photons adopt the phases of the a. c. modulating pulses, leading 
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to the effective controllable complex hopping constants in the rotating frame of Hs 1321133. 
Through the application of the developed three-tone PFC pulses to each of the grounding 
SQUIDs, every vertical hopping branch and every pair of horizontal hopping branches can 
be independently controlled by a modulating tone threaded in one of the grounding SQUIDs 
[35], implying that the artificial magnetic field for microwave photons with Landau gauge 


A = [0, Ay(x), 0], B = Be z = 


O’ 0, Q^Ay(x) 


( 6 ) 


can be created with in situ tunability. A further estimation demonstrates that the uniform 
hopping strength can be synthesized in the range |T/27 t| e [5,15] MHz [30, [3T, 35]. 

For the investigation convenience of the ESM physics, in what follows we endow a non¬ 
trivial ring geometry to the TLR lattice, i. e. an N x x N y square lattice with an n x x n y 
vacancy at its middle, as shown in Fig. [l](b) . Through the careful setting of the hopping 
phases, we penetrate a uniform effective magnetic flux 0 in each plaquette of the lattice and 
an extra a at the central vacancy. In Fig. [2|a) the energy spectrum of a finite lattice is cal¬ 
culated with N x = N y = 24 and n x = n y = 6. In the rational situation 0/27 t = p/q with p, q 
being co-prime integers, the unit-cell of the lattice is enlarged by q times, leading to q nearly 
flat magnetic bands and the fractal Hofstadter butterfly spectrum [3Z] (Fig. [2|a)). These q 
magnetic bands have nontrivial topological band structures, and between the q bands there 
exist ESMs traversing the q— 1 band gaps I31BLH3I- The lattice spectrums of the situations 
of interest p/q = 1/4 and p/q = 1/5 are shown in Figs. [2](b) and[2|c), respectively, where 
the flatness of the band steps and the stiffness of the connections between the steps imply 
the degeneracy of the Landau levels and the spectral location of the ESMs. 

Probing the ESMs: The spatial and spectral information. —Compared with 
fermionic electronic systems, the photonic nature of circuit QED allows multiple occupation 
of a particular mode at the same time and the non-equilibrium driving-dissipation competi¬ 
tion. Here we propose the following scheme of probing the topological nature of the ESMs. 
With the detailed modeling being discussed in Methods, we emphasize that the physics 
behind is that the exotic properties of the ESMs result in the novel steady states of the 
lattice, and the information of the ESMs can be extracted from the SSPNs of only few sites 
on the lattice versus the pumping frequency and the pumping sites. 

Firstly let us consider the single-site driving of a particular site r p described by 


n SP = Vsa rp e^ t + h.c., 


(7) 
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with Vs being the pumping strength and lisp being the detuning in the rotating frame of 
Hs- The SSPN on the pumping site rigp = (nj p a rp ) in the situation p/q = 1/4 and a = 0 
is numerically simulated based on equation ( [l6| ) and plotted in Fig. [3| In what follows we 
show that the spatial and spectral information of the ESMs can be distilled by measuring 
the dependence of the single-site SSPN (a£ a rp ) on Qgp and r p . 

If we choose r p = ro = (1,24) as an outer edge site (OES), significant ngp can be 
detected when Dgp falls in the 1st and 3rd gaps, indicated by the highlighted spectrum 
comb in Fig. |3](a) (for the even q = 4, the 2nd gap is closed as a Dirac point form). This 
can be attributed to the excitation of the outer ESMs. However, if Qgp is chosen deeply 
in the magnetic bands, n gp has bare value because in this situation Hsp can only excite 
bulk state modes (BSMs) which spread over the whole lattice, i. e. the weight of r p in 
the mode function becomes diluted. The situation of pumping an inner edge site (IES) 
r p = ri = (9,13) is similar, where the comb-like spectrum of ng P centralized in the band 
gaps can also be found in Fig. |3](b). Meanwhile, there are still several interesting differences. 
As the number of the IESs is smaller than that of the OESs, Fig. [3](b) contains fewer peaks 
than Fig. I») 0- In addition, different pumping strengths have been used in the numerical 
simulation of Figs. [3|a) and J^b) such that the obtained n r S p and n r s ' p are in the same region. 
This choice can also be traced back to the small number of the IESs which results in the 
concentration of the mode functions in the inner edge. Another observation is the opposite 
trends of Ugp and Ug P versus Qgp: In the 1st gap of Figs. [3]f a) and (b), the peaks of ngp/rigp 
increase/decrease with increasing Dgp. In contrast, when we set r p = rp = (5,13) as a bulk 
site (BS), the lattice will have detectable ngp iff Qgp falls in the magnetic bands. When 
we choose Osp in the band gaps, the lattice cannot be excited because no BSM spectrally 
populates in the band gaps and no ESM spatially populates in the bulk of the lattice (notice 
the marked window at the 1st and 3rd gaps in Fig. [3](c)) . 

The above illustration can be experimentally detected by the proposed measurement 
scheme sketched in Fig. [l](a) : A particular pumping site r p is capacitively connected to 
an external coil with input/output ports for pumping/measurement. The steady state of 
the lattice can be prepared by injecting microwave pulses through the input port for a 
sufficiently long time. During the steady-state period, energy will leak out of the r p th 
TLR from the coupling capacitance, which is proportional ay,, (''4,, °r p ) with the proportional 
constant being determined by the coupling capacitance. The target observable (aj p « rp ) can 
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therefore be measured by simply integrating the energy flowing to the output port in a 
given time duration. Actually, this measurement scheme has already been used in a recent 
experiment in which both the amplitude and the phase of a coherent state of a TLR were 
measured [32|. Here we emphasize that what we want to measure is the expectation value 
( a r p a r P ) ■ while the detailed probability of the multi-mode coherent steady state projected to 
the Fock basis is nevertheless not needed. It is this weak requirement that greatly simplify 
our measurement. 

We further calculate for each pumping situation a typical steady state photon distribution 
and display them in Figs. [3](d)— (f), respectively. In Figs. id) and |3^e) , the steady states 
correspond to the excitation of an inner or out ESM. The confinement and uniformity of the 
steady states clearly reflects the ESM mode functions localized and uniformly distributed on 
the edge, while the extended spatial distribution of the BS pumping steady state in Fig. [3|d) 
illustrates intuitively the difference between the BSMs and the ESMs. 

Probing the ESMs: Measuring the momentum. —While the single-site pumping 
provides a route of discovering the spatial and spectral properties of the ESMs, the more 
interesting physics comes from the multi-site inhomogeneous pumping, which proves to be 
an efficient method of measuring the momentums of the ESMs. We pump m consecutive 
OESs as 

m 

Wmp = Vm J2K ei{QMpt ~ jkp) + h - c 'L (8) 

3 = 1 

and investigate the summed SSPN %p = i ( ( 4, T,) on the m pumping sites. Here Vm 
is the homogeneous pumping strength, Ay is the phase gradient of the pumping between 
neighboring sites, Qmp is the frequency detuning in the rotating frame, and r ? for j = 
1, 2 ... m denotes the jth of the m pumping sites. 

Suppose f^MP matches the eigenfrequency of a particular outer ESM, there arises an 
interesting question that, how does timp depend on Ay? For a photon in that ESM, we can 
imagine its propagation around the edge with its ESM momentum k 0 (this can be verified 
by the discussion of the coherent dynamics in Discussion). Therefore, the reduced ESM 
mode function on the m pumping sites can be represented by a vector 

k 0 = (l,e iko ,e i2ko ...e i{m - 1)ko ) (9) 

where the e lko factor denotes the phase delay between two consecutive sites and the equal¬ 
weighting character of k 0 reflects the uniform spatial distribution of the ESM on the confined 
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edge (see Figs. |3](d) and[3](e)). It is this form of ko that inspires the inhomogenous multi-site 
pumping Hmp, which can be represented by another vector 

k P = (1, e ikp ,e 2ikp .. . e (m-1 )* fcp ). (10) 

Based on the above physical picture, we can conjecture that the maximum of ump will 
emerge at the point 

kp = ko, ( 11 ) 

where the excitations from the m pumping sites constructively interfere with each other. 
The dependence of «mp on kp and m is plotted in Fig. [4j where the positions of the peaks 
infer the value of ko- In addition, the full width of half maximum (FWHM) of the peaks 
decreases with the increase of m. This can be understood by considering the two extreme 
cases: If m = 1, there is certainly no peak because the steady state is independent on kp. 
Meanwhile, when all the OESs participate in the inhomogeneous pumping and the lattice 
size grows up, the inner product |kpk 0 | describing the interference between the pumping 
sites becomes nonzero iff kp = ko- In this situation the peaks in Fig. |4] approach a h-like 
function. As implied in Fig. [4| for a moderate m = 5 the FWHM is already sharp enough 
to discriminate the ESM momentums with a satisfactory resolution. 

Probing the ESMs: The integer topological invariants. —We further consider the 
measurement of the integer topological quantum numbers of the system. The topological 
property of a electronic Bloch band is captured by the quantized Hall conductance which 
turns out to be its Chern number [3j. This transport measurement is nevertheless inaccessible 
in circuit QED systems due to the absence of Fermi statistics. Meanwhile, the presence of 
the ESMs provides an alternative way of probing the topological invariants according to the 
bulk-edge correspondence [lj. In the rational situation 0/27 r = p/q, the eigenenergies of the 
ESMs can be represented by the zero points of the Bloch functions winding around the q — 1 
holes of a complex energy surface which correspond to the q — 1 gaps of the lattice |3]. The 
topological quantum numbers of the ESMs are given by the Diophantine equation 

h = s h q + t h p, \t h \ < q/2, (12) 

where h is the gap index and th , $h are integers. Especially, th is the topological winding 
number of the hth gap which is related to the Chern number (0 of the hth band as 


Ch — th ~ th-1- 


( 13 ) 


Measuring the winding numbers of the ESMs is thus equivalent to measuring the Chern 
numbers of the magnetic bands. For q = 4 we have C\ = C 3 = 1 and t\ = —t 3 = 1, while for 
q = 5 we have C 3 = —4 and Cj = 1 for the other four bands, and t\ = —£4 = 1 ,t 2 = —t 3 = 2 
(see Figs. 0b) andgc)). 

As the spatial configuration of the proposed lattice is equivalent to the Laughlin cylinder 
i (i- e. it can be regarded as the rolling of a two-dimensional simply-connected plane 
and the consequent threading with an effective magnetic flux a, see Fig. [ljc)) , the adiabatic 
pumping of the ESMs can be realized through the control of a in the central vacancy. Once 
a increases monotonically from 0 to 2ir, an integer number of ESMs will be transferred 
with the spectrum of the lattice returned to its original form. Such integer is exactly the 
winding number of the ESMs. In Figs. [5](a) and [5](d) the dependence of «mp on kp and a 
is numerically calculated for p/q = 1/4. Guided by the dashed white lines and the solid 
yellow arrows, the ESM peaks in the 1st and 3rd gaps move by one with the opposite moving 
directions, in agreement with the relation t\ = —t 3 = 1. This scheme is in principle general 
to measure integer topological invariants with any value: The peaks of the ESMs in the hth 
gap move by |G| during the variation of «, with the moving direction indicating the sign of 
th, and the Chern numbers C/, can be calculated from equation ( fT3[ ) after all £/, are obtained. 
The situation (f>/ 27r = 1/5 is also shown in Figs. [5^b) and[5](e), where the movements in the 
2nd and 3rd gaps cross two peaks, indicating f 2 = —£3 = 2. In addition, the opposite moving 
directions can also be observed in the pumping of the inner ESMs and the outer ESMs in 
the same gap (see Figs. and [5](f )). 


DISCUSSION 

Robustness against imperfection factors. —The imperfection in realistic experiments 
accompanies inescapably with the ideal scheme proposed above, including the diagonal and 
off-diagonal disorders u> r —> ui r + 5uj r ,T —>■ T + 47b r added into 74 t by the fabrication 
errors of the circuit and the low-frequency 1// background noises [14(3] , and the non-nearest- 
neighbor hoppings induced the residual long-range coupling between TLRs which are not 
presented in the proposed 74 t- Understanding these effects is thus crucial for our scheme. 
As briefly summarized below and detailedly studied in ref. [35], these imperfection factors 
lead to unwanted terms which are all much smaller than the band gaps of the lattice ( rT). 
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Moreover, some of the resulted effects can be further suppressed by the slight refinement of 
the proposed scheme. Such small fluctuations cannot destroy the topological properties of the 
ESMs because they are not strong enough to close and reopen the band gaps. Therefore, the 
presence of these imperfections can only renormalize negligibly the predicted results. From 
this point of view, our scheme enjoys the topological protection against the imperfection 
factors, which pinpoints its feasibility based on current level of technology. 

In particular, the fabrication error including the deviations of the realized circuit param¬ 
eters from the ideal settings (e. g. the lengths and the unit capacitances or inductances of 
the TLRs) leads to the disorder of the eigenmodes’ frequencies 5oj t , while the long-range 
coupling induced by the finite inductances of the grounding SQUIDs results in the next- 
nearest-neighbor hoppings through the current-division mechanism [38j 39J. As evaluated 
in ref. [35], these effects are both at the level of [10 -2 ,10 _1 ] T. In addition, these errors 
can be further corrected by simple revision of the proposed PFC scheme, including the re¬ 
fined choice of the TLRs’ lengths and the corresponding renormalization of the modulating 
frequencies of the grounding SQUIDs [35j. 

In actual experimental circuits, the 1// noise at low frequencies far exceeds the thermo¬ 
dynamic noise [10]. The 1// noise in superconducting quantum circuits can generally be 
traced back to the fluctuations of three degrees of freedom, namely, the charge, the flux, 
and the critical current. Due to its low frequency property, we can treat the 1 // noises 
as quasi-static, i. e. the noises do not vary during a experimental run, but vary between 
different runs. Firstly, the proposed circuit is insensitive to the charge noise as it consists of 
only TLRs which are linear elements and grounding SQUIDs which have very small charging 
energies and very large effective Josephson coupling energies. Such insensitivity roots in the 
same origin of the charge insensitivity of transmon qubits [41j and the flux insensitivity of 
the low decoherence flux qubit 021- Secondly, the flux 1// noises penetrated in the loops of 
the grounding SQUIDs shift the d. c. bias points of the grounding SQUIDs in a quasi-static 
way. The consequent effect is then the fluctuations 

Su r < 10" 3 T, 8Tr'r < 10~ 4 T, (14) 

where the detailed evaluation can be found in ref. [35]. These flux-noise-induced diagonal 
and off-diagonal fluctuations are both much smaller than the band gaps and the spectral 
spacing between the ESM peaks (~ 10 _1 T, see Fig. [ 5 ]). Such small fluctuations can thus 
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neither destroy the topological properties of the ESM nor mix the resolution of the ESMs in 
the SSPN measurement. Therefore, we come to the conclusion that our scheme can survive 
in the presence of the 1// flux noise. The influence of the critical current noise is similarly 
analyzed, with results indicating that the induced effects are even smaller than those of the 
flux 1// noises and can then be safely neglected [35, 43j. 

Coherent chiral photon flow dynamics. —A natural further step beyond the previous 
SSPN investigation is the study of coherent dynamics of the lattice, which is becoming 
experimentally possible due to the recent extension of the coherence times of superconducting 
circuits EM. Such investigation can offer an intuitive insight into the chiral property of the 
ESMs. We assume that the lattice is initially prepared in its ground state and then a driving 
%sp is added with r p being an edge state and ffgp being the eigenfrequency of an ESM. The 
time evolution of the lattice is calculated with screenshots of the photon flow dynamics 
shown in Fig. [6| The chirals of the ESMs result in the unidirectional photon flows around 
the edge with the directions determined by Dsp and r p . Chosen r p = (1,13) as an OES, 
the photon flow is clockwise/counterclockwise if ffgp falls in the lst/3rd gap (Figs. B&) and 
Be)), indicating the opposite chirals of the ESMs in different gaps. The energy separation 
of ESMs with different chirals can be understood by regarding the unidirectional flow as a 
rotating spin ra- Placed in an artificial magnetic field B, such spin has split energies with 
one spinning direction lower and the opposite direction higher. Moreover, it is observed 
from Figs. IN and [6](b) (see also Figs. B c ) and[6](d)) that the chiral of the inner ESMs is 
opposite to that of the outer ESMs in the same gap. This oppositeness can be explained 
by the spatial configuration of the lattice. As shown in the upper panel of Fig. B c )> by 
“tearing” the lattice apart we can get a simply-connected plane where there is only one edge 
existing. Now we consider the inverse, i. e. we “glue” the two sides marked by dashed 
lines together, and obtain the ring geometry shown in Figs. B a ) and the lower panel Of0c). 
During this gluing, the ESM flow (marked by the arrows) cancels itself on the glued sides, 
leaving two closed circulations with opposite directions. This “tearing-and-gluing” process 
can also be tested in our lattice configuration because the horizontal hopping branches can 
be adiabatically tuned on and off. Another interesting observation is the different velocities 
of the photon flows between Figs. Ba) and 13 c). Notice the pumping frequencies in these 
two subfigures are not mirror to each other, such difference in flowing velocity reflects the 
momentum difference between the corresponding ESMs. 
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The chiral flow in the presence of disorder and defect is also calculated. Here we assume 
that 8u r and T r ' r are normally distributed with a(5oj r ) = a(T r ' r ) = 0.057” much larger than 
those estimated in ref. [35 j. In addition, a 2 x 2 hindrance is placed on the upper outer edge 
with &U( 12-13,23—24)/T = 30. As displayed in Fig. [6j(e), the survival of the chiral photon flow 
under disorder and its circumvention around the hindrance clearly verify the topological 
robustness of the ESMs. 

We remark that the proposed unidirectional photon flow can also be detected through 
the photon-number measurement of only few sites neighboring to the pumping sites or 
the defect sites: We first establish the PFC process and pump the lattice with "Hsp for a 
duration less than the time scale during which the photon flow circulates around the whole 
edge loop and then remove them. The energy leaked out from the edge sites neighboring 
to the pumping sites can be observed, which is proportional to the photon numbers stored 
in those TLRs. From Fig. [6] we expect that the injected photons tends to flow towards 
the sites on a particular direction with phase delay while leaves the sites on the opposite 
direction negligibly excited [21(331. The comparison of the measured photon numbers on 
the opposite directions thus reveals the chiral property of the ESMs. The circumvention of 
the photon flow around the hindrance can be detected in a similar manner by measuring 
the sites neighboring to the hindrance sites. 

Extension in the future. —The flexibility of the proposed PFC method is not limited 
by the square lattice configuration in this paper. For instance, a brick wall lattice (i. e. 
a stretched honeycomb lattice) can be straightforwardly obtained by closing some of the 
vertical hopping branches in the original square lattice. This generalization may pave an 
alternative way to the study of photonic graphene j33j. In addition, while the cross-talk 
between the diagonal next-nearest-neighbor TLRs is suppressed in our scheme due to fre¬ 
quency mismatch, it can indeed be opened by adding another tone to the modulating pulses 
of the SQUIDs. This may offer potential facilities in the future study of anomalous quantum 
Hall effect in the checkerboard lattice configuration ESI- 

In recent research, lattice configurations supporting a dispersionless flat band have been 
investigated extensively, including the Lieb and the Kagome lattices [25] . These band struc¬ 
tures provide an idea platform of achieving strongly-correlated phases as the kinetic energy 
is quenched 051116]. These lattice configurations can also be synthesized by the variation 
of the proposed square lattice. While the photonic topological insulator considered in this 
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paper can be understood in the single-particle picture, the introduction of interaction signif¬ 
icantly complicates the problem and may lead to much richer but less explored physics. On 
the other hand, with the demonstrated strong coupling between TLRs and superconducting 
qubits [7, S] (and also atomic system, see ref. [42]), the Bose-Hubbard [Tbl fl8] and Jaynes- 
Cummings-Hubbard nonlinearities na H9] can be incorporated into the proposed lattice. 
A further research direction should therefore be the implementation of photonic fractional 
Chern insulators and the understanding of strong correlation in the proposed architecture 
and its potential hybrid-system generalizations which may utilize the advantages of different 
physical systems mm- 

CONCLUSION 

In conclusion, we have proposed a method of implementing topological photonics in a 
circuit QED lattice. The effective magnetic field for microwave photons can be synthesized 
through the proposed parametric approach, and the topological properties of the ESMs can 
be extracted from the steady states of the lattice under pumping. Moreover, being flexible to 
incorporate effective photon-photon interaction, our proposal may offer a new route towards 
the investigation of nonequilibrium photonic quantum Hall fluids in on-chip superconducting 
quantum circuits. Taking the advantage of simplicity in setup and topological robustness 
against potential imperfections, the realization of this scheme is envisaged in the future 
experiments. 


METHODS 


Steady state of the lattice. —We explicitly consider the pumping of the lattice de¬ 
scribed by pfae - * m + h.c. in the rotating frame of "Hs • Here P and a are the vectors 
composed of the pumping strengths and the annihilation operators of the lattice sites, re¬ 
spectively, and fl is the monochromatic detuning playing the role of Fermi surface. In the 
presence of dissipation, the evolution of the lattice is described by the master equation 


dp 
d t 


[a- f2Z)a + P*a + aJP, p] 


+ 


n r (2 a r pa\ — a\a r p — paja r ) , 


(15) 
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where p is the density matrix of the lattice, n r is the decay rate of the rth TLR, and the 
matrix B is defined by a^Ba = 'Hj. As a linear system (i. e. there is no photon-photon 
interaction), the lattice can be described in the picture of multi-mode coherent state, and 
its steady state can thus be determined by 


,d(a) 
d t 


B -9,1 - -iJC 
2 


(a) + P = 0, 


(16) 


where /C is the diagonal matrix of the TLRs’ decay rates. From equation (16) our idea 


emerges that, 9 can be used to select the mode we are interested in, and the information of 
that mode can be extracted from the dependence of (a) on P. In addition, as the pumping 
sites are coupled to external coil, they suffer more severe decoherence than the other conven¬ 
tional sites. Therefore, we set the decay rates of the pumping sites as uniformly /v rp /2 tt = 2 
MHz and those of the conventional sites as uniformly K r /27r = 100 kHz throughout the 
numerical simulation of this paper, i. e. there is a 20 times difference between them 
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FIG. 1. (a) Sketch of the square TLR lattice with the four colors (yellow, blue, green, and red) 

denoting the different lengths of the TLRs and the black dots representing the grounding SQUIDs 
(see the lumped circuit magnified representation at the right side). Each TLR plays the role of a 
photonic site (the large colored rounds) and the effective hopping between them (the dotted-dash 
lines) can be induced through the dynamic modulation of the SQUIDs. The pumping and the 
consequent steady-state measurement can be performed through the external coil connected to 
the pumping site(s) (lower right), (b) Configuration of the proposed lattice. The colored rounds 
and the solid lines label the TLRs with corresponding lengths and the photon hopping branches, 
respectively, (c) Spatial geometry of the lattice. The lattice shown in (b) can be obtained from 
the gluing of a simply-connected plane by the two dashed sides (the upper panel). Through this 
process the opposite chirals of the inner and the outer ESMs (the arrows) are formed (the lower 
panel). 
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FIG. 2. (a) Hofstadter butterfly spectrum of the proposed lattice with N x x N y — 24 x 24, 

n x x n y — 6 x 6, and a/2n — 0, where energy is in units of T. The situations of the rational 
effective magnetic fields cj)/2n — 1/4 and </>/ 2n — 1/5 are denoted by the red and blue lines, 
respectively. Their eigenenergies are shown in (b) and (c), with m the index labeling the 540 
eigenvalues from smallest to largest. The band gaps are highlighted with their topological winding 
numbers marked. 
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FIG. 3. Steady state of the proposed TLR lattice with </>/27r = 1/4, a = 0, and T /27r = 10 MHz. 
The SSPN on the pumping site rigp versus flsp are displayed in (a), (b), and (c), with r p = [(1, 24), 
(9,13), (5,13)] and Vs/T = [0.5, 0.25, 0.23], respectively. The representative SSPN distributions 
on the whole lattice are presented in (d)—(f), with the pumping frequencies fisp/7~ = [1.47, 1.97, 
2.69] marked by the corresponding red lines in (a)—(c). 
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FIG. 4. nMP versus kp and m with Vm/T = 0.1 and Qmp/T = 1.67. The pumping sites start 
from the OES site ([7 — m/2 \, 24) and end at ([6 + m/2j, 24). The other parameters are the same 
as those in Fig. [3j 
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FIG. 5. Adiabatic pumping of the proposed lattice represented by timp versus (Imp and fcp. Here 
we set (f)/2n = 1/4 for (a), (c), (d), and (f), and (f)/2n = 1/5 for (b) and (e). The panels from 
top to bottom in each of the subfigures correspond to a/2n = 0, 1/4, 1/2, 3/4 and 1, respectively. 
For (a), (b), (d), and (e) the OESs (4,1)—(8,1) are pumped, while for (c) and (f) the IESs 
(9,9)—(13,9) are pumped. Notice the ranges of the pumping frequency are chosen in a mirror 
form for the upper and lower subfigures. The other parameters are the same as those in Figs. [3] 
and|U 
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FIG. 6. Chiral flow dynamics in the presence of single-site pumping. We set r = (1,13), Vs/T — 2 
for (a), (c), and (e), and r = (9,13), Vs/T — 1.5 for (b) and (d). The synthesized magnetic field is 
set as <fi/2n = 1/4 and a/27r = 0. In addition, Qsp /T is chosen as [—1.76, —1.97, 1.47, 1.97, —1.75] 
for (a)-(e), respectively. The times of the panels from left to right are arithmetic progressions 
with the first terms Ti/2n = [6, 1, 6, 1, 15]T -1 and the common differences AT/2tt = [13, 5, 13, 5, 
2.5]T -1 for (a)-(e), respectively. Especially, in the calculation of (e) the effect of lattice disorder 
and defect is incorporated (see the main text). The other parameters are set the same as those 
used in Fig. [3j 
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